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The hadronization scheme for parton transport in relativistic heavy ion collisions is considered 
in detail. It is pointed out that the traditional scheme for particles being freezed out one by 
one leads to serious problem on unreasonable long lifetime for partons. A super-cooling of the 
parton system followed by a collective phase transition is implemented in a simple way. It turns 
out that the modified model with a global phase transition is able to reproduce the experimental 
longitudinal distributions of final state particles better than the original one does. The encouraging 
results indicate that a relevant parton transport model for relativistic heavy ion collision should 
take equilibrium phase transition into proper account. 
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I. INTRODUCTION 

The quark-gluon plasma (QGP) is expected to be 
formed in heavy-ion collisions at the Relativistic Heavy- 
Ion Collider (RHIC). So far very interesting experimental 
data have been collected There are strong evidences 
on the deconfinement of QCD vacuum and the appear- 
ance of a (locally) thermalized partonic system at the 
early stage of collision. Both theoretical and experimen- 
tal investigations are stimulated on the evolution of the 
partonic system and the way it transforms to final state 
particles. 




FIG. 1: The different stages of heavy-ion collisions. 



The reaction process of relativistic heavy ion collision 
can be divided into several stages (see figure At first 
the primary interaction creates many partons. These 
partons are certainly far from equilibrium and cascade 
through interaction among each other. If the collid- 
ing nuclei is heavy enough and the incident energy is 
high enough the partonic system is expected to reach 
equilibrium or local equilibrium, forming quark-gluon 
plasma. The system then expands and the temperature 
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decreases, passing via a phase transition or crossover into 
the hadron gas stage. Hadrons continue to cascade until 
they freeze out from the collision region, forming final 
state particles. 

At present there is no unique theory that can describe 
the reaction process as a whole. Different theoretical ap- 
proaches are applied to different stages [2|. The primary 
interaction, creating many partons is often described by 
eikonalized parton model Jjl, Gribov-Regge theory Q or 
parton saturation model |5j- Parton and/or hadron cas- 
cade is referred to as the solution of non-equilibrium 
transport equation R| and is usually realized by Monte 
Carlo models B The hydrodynamics 01 or ther- 
mal models [Tl| are employed to deal with the parton 
and/or hadron system in local or global equilibrium. 

Usually, the hydrodynamics or thermal models take 
(local) equilibrium as model assumption and do not an- 
swer the question on how the system arrive at (local) 
equilibrium. The global properties of the system are 
the main issues considered in the model, while the de- 
tailed evolution of the constituents — — partons and/or 
hadrons is not taken into account. 

On the contrary, the transport models follow the evo- 
lution of the partons and/or hadrons in detail through 
considering their interaction — — elastic and/or inelas- 
tic scattering, or cascade. Usually, the model allows the 
cascade continues to go on until the interaction ceases 
and then the parton or hadron freeze out from the sys- 
tem. A noticeable common property of such model is 
that, freeze out of parton or hadron occurs particle by 
particle . Each particle has its own freeze-out time. Such 
an approach is acceptable for a hadronic system, where 
hadrons after ceasing to interact with other hadrons fly 
away from the system freely toward the detectors, but 
will cause serious problem for a partonic system. 

Partonic and hadronic systems are of different phases. 
In general, a phase is a portion of a system that is uni- 
form and has a definite boundary |12|. In particular, the 
partonic and hadronic phases exist in different vacua — 
- the former is in perturbation QCD vacuum and the 
latter in physical vacuum. Therefore, the transition be- 
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tween partonic and hadronic phases should be a collec- 
tive, thermal-equilibrium phenomenon, accompanied by 
a vacuum transform, which could not be realized in a 
particle-by-particle way. 

In the transport models presently on market, the 
hadronization is realized parton-wise instead of collec- 
tively, and therefore, in these models there is only 
hadronization but no partonic to hadronic phase transi- 
tion in the strict sense. To let a few partons live for a very 
long time is inconsistent with the general believe that the 
existing time of QGP is about 1-5 fm/c and hadron 
freeze out at about 20-40 fm/c 01- Even more seriously, 
when most of the partons have already hadronized, the 
system is dominated by hadrons and the corresponding 
vacuum is a physical one instead of a QCD perturbation 
one. To let some partons survive and fly freely in such a 
circumstance is highly unphysical. 

The aim of the present paper is to discuss this prob- 
lem in detail. We will take as example a presently avail- 
able model, that has parton transport implemented. The 
temperatures of the system at different time will be ex- 
tracted using thermal-equilibrium transverse mass distri- 
bution. The partonic system will be allowed to somewhat 
super- cooled, i.e. the temperature is allowed to decrease 
to lower than the expected phase-transition temperature. 
At a certain point all the partons remaining in the sys- 
tem are forced to coalesce, forming hadrons. Then the 
hadrons start to cascade toward freeze out. Thus the dif- 
ficulty of long-life parton is overcome in a simple manner. 
The phenomenological consequences of such an approach 
will be presented and compared with existing experimen- 
tal data. The possible reason for the improvement of the 
present approach with global hadronization in compari- 
son with the original model with long-life partons will be 
discussed. 

The layout of the paper is as the following. A short 
introduction to the AMPT model, used as an example 
for models with parton transport is given in Section II. 
A detailed analysis about the parton and hadron time 
evolution in AMPT is then presented in Section III. A 
super-cooling of the parton system followed by a collec- 
tive hadronization scheme is proposed in Section IV to- 
gether with the phenomenological consequences on the 
final hadron distribution and elliptic flow. Section V is 
conclusion and discussion. 



then modeled by the ZPC parton cascade model. At 
present this model includes only parton-parton elastic 
scattering with cross section 



where the screening mass fi is taken to be an input pa- 
rameter of the model to obtain the desired total cross 
section. Two partons will undergo scattering when the 
closest distance between them is smaller than y^a/n. 

There are two versions of AMPT model. In the de- 
fault AMPT, after ceasing interactions minijet partons 
are combined with their parent strings to form excited 
strings, which are then converted to hadrons according 
to the Lund string fragmentation model. While in the 
AMPT with string melting, the strings in the initial con- 
ditions are melt to partons first and then interactions 
among all the partons are again modeled by ZPC. Af- 
ter partons stop interacting, a simple quark coalescence 
model is used to combine the two nearest partons into 
a meson and three nearest quarks (antiquarks) into a 
baryon (antibaryon). Scatterings among the resulting 
hadrons are described by a relativistic transport (ART) 
model [Hi| which includes baryon-baryon, baryon-meson 
and meson-meson elastic and inelastic scatterings. 

It turns out that the default AMPT (vl.ll) is able to 
give a reasonable description on hadron rapidity distri- 
butions and transverse momentum spectra observed in 
heavy ion collisions at both SPS and RHIC. However, 
it fails to reproduce the experimental data about ellip- 
tic flow and two-pion correlation function. On the other 
hand, the AMPT model with string melting (v2.11) can 
well describe the elliptic flow and two-pion correlation 
function [13 H3 but agree baddly with the hadron ra- 
pidity and transverse momentum spectra. 

In the following we will utilize the AMPT with string 
melting — — AMPT v2.11 to generate Au-Au central col- 
lision events at ^/snn = 200 GeV. The impact parameter 
is in the range b < 3 fm and the parton cross section is 
taken to be 10 mb. 



III. THE TIME EVOLUTION OF PARTONS 
AND HADRONS IN AMPT 



II. A BRIEF INTRODUCTION TO AMPT 

The AMPT model 9] is based on non-equilibrium 
transport dynamics. It contains four main components: 
the initial conditions, partonic interactions, conversion 
from the partonic to the hadronic matter and hadronic 
interactions. The initial conditions, which includes the 
spacial and momentum distributions of minijet partons 
from hard processes and strings from soft processes, are 
obtained from the HIJING model in which eikonized par- 
ton model is employed. The time evolution of partons is 



In AMPT v2.11, each initial parton has a formation 
time given by tf = Eh /m| h with Eh, tut^h the energy 
and transverse mass of its parent hadrons. After this 
formation time the partons start to scatter with each 
other and when a parton no longer scatters with any 
other parton, it will hadronize. Thus each parton has its 
own hadronization time. 

In Fig. [21 are shown the percentages of parton and 
hadron, respectively, at different time after the collision. 
It can be seen from the figure that a few hadrons (about 
4%) have already emerged at t < 5 fm/c. At this time, 
partons dominate, and the system as a whole is in the 
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deconfined phase, located in the perturbation QCD vac- 
uum, with a few hadrons vaporized out. 

As time increases, the number of partons decreases 
while that of hadrons increases. In this process a par- 
ton transforms to hadron when and only when it ceases 
to interact with other partons. Thus a part of parton 
survive up to an unreasonable long time, e.g. t ~ 100 
fm/c. This is the difficulty common for this kind of mod- 
els mentioned in the Introduction. 
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fm/c, slightly above the predicted critical temperature 
T c = 170 MeV 113, while at t = 5 fm/c the temper- 
ature has already arrived at about 110 MeV, i.e. in the 
first few fm/c the temperature decreases rapidly and soon 
becomes lower than the expected phase transition tem- 
perature. 
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FIG. 3: The transverse mass distribution (histograms) for 
d-quarks at two different time t — 0.5 and 5 fm/c . The lines 
are the fit to Eq. (3) 



FIG. 2: The percentage of partons and hadrons, respectively, 
in AMPT v2.11 for v /s^zv = 200 GeV Au-Au central collisions 
with b < 3 f m and parton cross section 10 mb. 



For simplicity, we assume that the system arrives at 
equilibrium globally instead of locally and the equilib- 
rium is examined at different time t instead of at different 
intrinsic time r. Under this assumption the temperature 
of the system at different time is determined by thermo- 
dynamics. 

Suppose the expansion of system is a slow and quasi- 
static process, then at each time step the system can be 
regarded as an equilibrium thermal system with temper- 
ature T. The invariant momentum distribution of parti- 
cles emitted from the thermal system is |19| 
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d 3 p dym^drri'Tdij) (27r) ; 



Ee - ( E- M )/T^ (2) 



where g is the particle spin-isospin degeneracy factor, fi is 
the chemical potential, V is the system volume, E is the 
energy of the particle, y is rapidity and wt = \J m 2 + p\ 
is the transverse mass. Integrating with respect to y and 
<fr, we get the transverse mass distribution 
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in which K\ is the hyperbolic Bessel function of rank one. 
By fitting the transverse mass distribution, the temper- 
ature of the system can be extracted. As example, in 
Fig's. |3| are shown the transverse mass distributions for 
d quarks at two different time t = 0.5 and 5 fm/c and 
the fit to Eq. (3). In the present work we assume the 
system to be static, and omit the effect of radial flow on 
the fitted temperature. 

It can be seen from Fig's. El that the temperature 
of the parton system is about 180 MeV at t = 0.5 
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FIG. 4: The temperature evolution extracted from parton 
and hadron transverse mass spectra in AMPT v2.11 with 
a collective phase transition implemented at t — 5 fm/c for 
^/sjvjv = 200 GeV Au-Au central collisions with b < 3 fm and 
parton cross section 10 mb. 



IV. A HADRONIZATION SCHEME WITH A 
SUPER-COOLING FOLLOWED BY A 
COLLECTIVE HADRONIZATION 

We regard the temperature of parton system decreas- 
ing rapidly to lower than the critical temperature as a 
supercooling effect, i.e. after the formation of QGP, the 
system lowers its temperature by expanding and evapo- 
rating hadrons, arriving at a temperature lower than the 
phase transition temperature. Then at a certain point all 
the left partons are coalesced to hadrons. The duration 
of the supercooling state is taken as a model parame- 
ter and in the present work we take t = 5 fm/c (the 
corresponding parton temperature is about 110 MeV) to 
be the beginning of the phase transition. At this time, 
all the remaining partons stop to interact and start to 
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hadronize. The temperature of the resulting hadron sys- 
tem is higher than that of the parton system due to the 
release of latent heat, and then the temperature decreases 
again through expansion, cf. FigQ] 

Thus we have implemented an equilibrium phase tran- 
sition with supercooling to the transport model in a very 
simple way. Our purpose is to see how the phase transi- 
tion affect the final state hadron distributions. 

The rapidity distribution is 



dN 
dy 



1 An 
Ay' 



(4) 



where 7V 0V is the number of events, Ay is the width of 
rapidity bin, An is the number of particles inside the ra- 
pidity bin. In Fig's. |3]are shown the rapidity distributions 
for charged particles, pions, kaons, net-protons, protons 
and anti-protons. The solid lines represent the results 
of AMPT v2.11 with parton- wise hadronization and the 
dotted lines are that with collective phase transition im- 
plemented in the above-mentioned way. Full circles are 
PHOBOS 0-6% centrality data and BRAHMS 0-5% cen- 
trality data. The data of dN^/drj are with both statisti- 
cal and systematic errors and the other data are with only 
statistical ones. The figures show that the rapidity distri- 
butions of the model with phase transition implemented 
have a better description of the experimental data, espe- 
cially on the kaon and proton rapidity distributions. As 
about the transverse distributions we found that they are 
almost unaffected by the implementation of global phase 
transition 

The success of AMPT with string melting is that it is 
able to describe the elliptic flow data very well through 
adjusting the parton cross section. So, whether the im- 
plementation of a global phase transition will destroy this 
agreement between model and data is a natural question. 

In Fig. are shown the elliptic flows V2 as a function 
of the number of participants calculated from the AMPT 
(with string melting) model with parton-wise hadroniza- 
tion and with collective phase transition, respectively. 
The parton cross sections are taken to be 3 mb and 10 
mb. It is clear from the figure that the V2 with collective 
phase transition implemented has almost the same shape 
as that with parton-wise hadronization. The implemen- 
tation of global phase transition affects the elliptic flow 
little and preserves the agreement between model and 
experimental data. 



CONCLUSION AND DISCUSSION 




FIG. 5: The rapidity distribution for ^/s^v = 200 GeV 
Au-Au central collisions. The solid lines are AMPT v2.11 
with parton-wise hadronization and the dashed lines are that 
with collective phase transition implemented. The impact 
parameter is b < 3 and parton cross section 10 mb. The 
dots are data from PHOBOS 6% and BRAHMS 5% central 
collisions. 



decreasing of temperature the parton phase is allowed to 
hadronize as a whole after a supercooling stage. It turns 
out that the modified model with a global phase tran- 
sition inherits the success of the original one in elliptic 
flow and is able to reproduce the experimental longitu- 
dinal distributions of final state particles better than the 
original one does. 
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FIG. 6: Centrality dependence of charged hadron elliptic flow 
for y'sjvjv = 200 GeV Au-Au central collisions. The solid 
lines are the AMPT v2.11 with parton-wise hadronization 
and dashed lines are that with global phase transition imple- 
mented. The dots are data from PHOBOS experiments |2^l . 
The impact parameter b < 3 and parton cross section is lOmb 
and 3mb, respectively. 



In the traditional transport models the particles in the 
system are freezed out one by one, each particle has its 
own freeze-out time. Extending such an approach to par- 
ton transport leads to serious problem on unreasonable 
long lifetime for partons. 

In order to avoid this problem, we assume the system 
to have reached global equilibrium and extract tempera- 
ture from the system by thermodynamics formula. As the 



In order to see why the model with collective phase 
transition can describe the experimental data better, the 
rapidity distributions of partons right before hadroniza- 
tion are plotted in Fig \7\ Comparing with the model 
with parton-wise hadronization, the parton transport in 
the model with collective phase transition is truncated, so 
there are fewer partons in the mid-rapidity region and the 
distribution peaks at regions with large absolute values 
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of rapidity. This effect results in the hadron distribution 
in mid-rapidity being lower than that from the original 
model, and thus approaching the experimental data. 

The elliptic flows in parton cascade models are built-up 
very early [ljj [23| (less than 5 fm/c), thus the truncation 
of parton transport at 5 fm/c in the present work does 
not affect the elliptic flow. 



FIG. 7: The rapidity distribution of partons right before 
hadronization in AMPT v2.11 with parton-wise hadronization 
(solid line) and with global phase transition (dashed line.) 
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We have proposed a super-cooling followed by a global 
hadronization as a prototype of the thermal-equilibrium 
phase transition from parton transport to hadron gas. 
Though our method for the implementation of phase 
transition seems to be very crude comparing to the real 
process in relativistic heavy ion collisions, the encour- 
aging results indicate that a relevant parton transport 
model for relativistic heavy ion collision should take equi- 
librium phase transition into proper account. 
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